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Abstract 

Inelastic neutron scattering (INS) experiments have been performed on the CU3 triangular molec- 
ular nanomagnet using powder samples. In the medium resolution INS spectrum measured, there 
are two peaks at tkij = 0.5 and 0.6 meV. Comparing the observed Q dependences of these peaks 
with calculations, these two INS peaks originate from the CU3 cluster. From the observed peak 
position, width, and intensity, we have determined the optimum parameters of the spin Hamil- 
tonian consisted of the antiferromagnetic exchange and the Dzyaloshinsky-Moriya interactions, 
which can also reproduce the magnetic susceptibility measurement. In addition, we have directly 
observed that the ground state quartet is split into two doublets with the energy separation of 
0.103 meV using high-resolution neutron spectroscopy, which exactly corresponds to that expected 
from the optimum parameters obtained from the medium resolution experiment. The temperature 
dependences of the integrated intensities of the 0.5 and 0.6 meV peaks are well reproduced by the 
Boltzmann distribution of the energy levels of the model Hamiltonian below 10 K. Furthermore, the 
inelastic peaks were visible even at very high temperatures as 50 K. This indicates extraordinary 
weak coupling between phonons (or any other perturbations) and spin states in the CU3 cluster, 
compared to the other known molecular nanomagnets. 
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I. INTRODUCTION 



Quantum phenomena in magnetism sometimes appear in the macroscopic measurement. 
Especially, molecular nanomagnets jlj] have been providing rich playgrounds to investigate 
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new types of such quantum phenomena. For instance, Mn 12 2] and Fe§ [3|, which are 
the most widely investigated spin clusters realized as the single molecular magnet, show 
a superparamagneticlike behavior and a quantum tunneling between the total spin states. 
Cr 7 Ni [4], where one Cr ion (S = 3/2) of the Cr 8 -ring [5| is substituted by a Ni ion (S = 1), 
is also a good example, showing a quantum coherence between the total spin states at the 
level- crossing field. 

Another interesting issue is a magnetization reversal with a pronounced hysteresis ob- 
served in V15 p). This phenomenon can be quantitatively explained by the dissipative 
two-level model consisted of a Landau-Zener-Stiickelberg (LZS) transition and a phonon 
bottleneck effect. For applying the LZS transition, the Dzyaloshinsky-Moriya (DM) interac- 
tion J7|, |8] (or the hyperfine interaction [9]) performs an important role to realize the avoided 
level crossing of the S to tai = 1/2 ground state doublets at zero field, where S to tid represents 
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FIG. 1: (a) A structure of the CU3 spin cluster, where each circle represents the Cu 2+ ion. The 
distances between Cu 2+ ions di2, c?2,3, and ds t i are written in the text, (b) Schematic view of the 
energy levels of the model Hamiltonian obtained by the optimum parameters given in Eq. ©; all 
lines stand for each doublets and the number next to the lines is the label of the energy levels. Hujo 
represents the splitting of the Stotai = V^ quartet by the DM interaction, whereas tvjj\ and huj2 
correspond to the INS peaks at 0.5 and 0.6 meV in the INS spectra measured at HER as shown 
in Figs. [5] and O The values ahead of arrows are the ratio of the calculated INS intensities of each 
excitations at \Q\ = 0.95 A -1 divided by that of hwo. 



the total spin of the cluster in this paper. Indeed, the zero field splitting of the ground state, 
which is a direct consequence of the DM interaction, was observed by the inelastic neutron 
scattering (INS) measurement ld j. In addition, from the result of changing the thermal 
coupling to the reservoir in the magnetization measurement, the phonon bottleneck effect is 
revealed to be the origin of the hysteresis [ll| . NaFeg 12] is another example in this cate- 
gory, which shows hysteresis loop near the level crossing field in the magnetic torque curve 
caused by the phonon bottleneck effect. On the other hand, the long life time of the spin 
state is the origin of the hysteresis in the V6 cluster 13[, which also shows the magnetization 
reversal with hysteresis. 

Recently, it has been found that the spin trimer clusters such as V 3 [14| . Cu 3 As 15| . 



f-step 



15| is 



and CusSb [16|, [17J show a similar but more intriguing phenomenon, which is a nal 
magnetization change with a milli-second order hysteresis. The DM interaction 7J, 
again the key to realize the avoided level crossing of the total spin, when the ground state 
changes to the Stotai = 3/2 state from the Stotai = 1/2 state. The long life time of the spin 
state is inferred to cause the hysteresis in the triangular cluster system. To date, there are 

nn 

only few reports [15|, H7J on the direct observation of the DM interaction, although we think 
the data are not sufficiently convincing. In addition, none for the origin of the long life time 
for these spin trimer system. These two facts are necessary to confirm the validity of the 
model for the half-step change with hysteresis in the pulsed-field magnetization curve. 



Nai2[Cu3(SbWg 033)2 (H 2 0)3]-46H 2 (hereafter, CU3 in short) 18] is one of such spin 
trimer clusters. Several measurements such as magnetic susceptibility, magnetization 
in pulsed field, electric spin resonance, and nuclear magnetic resonance have been per- 
formed [la, ll7J. Three Cu 2+ ions are placed at the distances of d\p = d 2 $ = 4.871 and 
t/31 = 4.772 A as shown in Fig. [Ha). Cu 2+ ions with 5* = 1/2 are coupled by the antifer- 
romagnetic exchange interaction, as indicated by the magnetic susceptibility measurement. 
The exchange path is thought to be Cu-O-W-O-W-O-Cu 16] • There is no inversion sym- 
metry at the center of any two Cu 2+ ions, suggesting the existence of the DM interaction in 
addition to the super exchange interactions. Consequently, the following spin Hamiltonian 



has been proposed [171 ]: 



3 x,y,z 



3 

+ HB^Si-g- B, (1) 



8=1 



where Jf i+1 and Df i+1 are the a-component of the exchange and DM interactions between 
the i-th and (i + l)-th Cu 2+ ions, and /xb is the Bohr magneton. The reported energy 
level scheme consists of the ground state with S tota [ = 1/2, the first excited state with 
Stotai = 1/2 at 100 /ieV hi gher than the ground state, and two weakly split S to tai = 3/2 



doublets at 580 /ieV higher [17 ] 



The splitting of the S to tai = 1/2 quartet into two doublets is due to the DM interac- 
tion ]19]. Thus, the DM interaction can be observed directly by detecting the excitation at 
hujQ. On the other hand, origin of the long life time of the spin state may be elucidated by 
directly observing the life time of the excitation levels, i.e., the broadening of the excitation 
peaks in INS spectra. In this paper, we will discuss these features using neutron scattering 
spectroscopy. First, we obtain the parameters of the model Hamiltonian of the CU3 spin 
cluster by observing the excitations to S tota i = 3/2 from S tota i = 1/2. Then, we observe 
the DM interaction microscopically by detecting the splitting of the Stotai = 1/2 quartet. 
Finally, we investigate the temperature dependence of the INS peaks to get an insight into 
the origin of the long life time of the spin state in CU3. 



II. EXPERIMENTAL DETAILS 

The powder sample was prepared by the procedure reported in the Refs. [l8| and 
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, and 



the deuterated powder sample was also prepared for the high-energy-resolution INS mea- 
surement. Magnetic susceptibility measurement of 45.3 mg non-deuterated powder sample 
was performed with a SQUID magnetometer in the temperature range of 1.8 < T < 300 K. 
A part of INS experiments with about 18.2 g non-deuterated powder sample was per- 
formed using the triple-axis spectrometer ISSP-HER, installed at the JRR-3 research reactor 
(Tokai, Japan). We have employed vertically focusing monochromator to select incident neu- 
tron wavelength, whereas double focusing (i.e. both horizontal and vertical focusing) tech- 
nique was used for the analyzer. Pyrolytic graphite (PG) 002 reflections were used both for 



the monochromator and analyzer. The spectrometer was operated in the fixed-final-energy 
mode with Ef = 2.4 meV, resulting in the instrumental resolution of 61 /xeV (FWHM, or full 
width at half maximum) at the elastic position. The resolutions at Huo = 0.49 and 0.60 meV 
are estimated as 68 and 71 /xeV (FWHM), assuming the Cooper-Nathans type resolution 
function [21] . The higher-energy harmonics were eliminated using the cooled Be filter. The 
non-deuterated powder sample was sealed in the aluminum sample can filled with the 4 He 
exchange gas, which is set to a closed-cycle 3 He refrigerator with the lowest working tem- 
perature of about 0.7 K. An ILL Orange cryostat was also used in the measurement of Q 
dependences as shown in Figs. 11(b) andH^d). 

An INS experiment was also performed using the disk chopper time-of-flight spec- 
trometer DCS installed at NIST Center for Neutron Research (Gaithersburg, USA) with 
E\ = 1.0 meV. The resolution at elastic position was obtained as 18.7 /xeV (FWHM), and 
the resolutions at hw = —0.1 and 0.1 meV were estimated as 22.6 and 15.1 /^eV (FWHM), 
respectively [22j. The deuterated powder sample of about 4.7 g was put in the aluminum 
sample can, and set to the ILL Orange cryostat, with which the lowest working temperature 
is 1.5 K. 

III. EXPERIMENTAL RESULTS 

An INS spectrum at \Q\ = 0.95 A" 1 and T = 0.71 K measured at HER is shown in 
Fig. |2fa). Although the strong incoherent scattering from hydrogen centered at the elastic 
position was observed as the background, there is an obvious INS peak around ftw = 0.6 meV. 
To focus on the INS peak, we plotted INS spectra at \Q\ = 0.95 A -1 , and T = 0.71 and 
2.62 K in the range of 0.3 < Hoo < 0.8 meV in Fig. Hl^b). The spectrum at 0.71 K exhibits a 
weak peak at 0.5 meV as well as the pronounced peak at 0.6 meV. At 2.62 K, the intensity at 
0.5 meV increases, whereas the intensity of the 0.6 meV peak decreases. The INS spectra of 
\Q\ = 0.95 A -1 at various temperatures measured at HER are also shown in Figs. EJ^c) — [2](j ) . 
The intensity of the 0.5 meV peak increases as temperature increases up to 2.62 K and 
decreases as temperature increases above 2.62 K, whereas the intensity of the 0.6 meV peak 
monotonically decreases as temperature increases. As seen in Fig. |2} the magnetic signals 
are rather smeared by the strong incoherent background. Therefore, to see the magnetic 
component clearly, the incoherent intensity was removed by the following procedure. The 
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FIG. 2: (a) The INS spectrum at \Q\ = 0.95 A" 1 and T = 0.71 K measured at HER is shown. 
Note the log scale for the vertical axis. A solid line is the fitting result of Eqs. (J2I) , (J3j) , and fl3J). 
(b) The INS spectra of \Q\ = 0.95 A" 1 at T = 0.71 and 2.62 K. The error bars of the observed 
data in both figures are within each marks, and the error bars in all figures represent one standard 
deviation. The INS spectra of \Q\ = 0.95 A" 1 at (c) T = 0.71 K, (d) 1.40 K, (e) 2.62 K, (f) 4.03 K, 
(g) 6.05 K, (h) 10.17 K, (i) 30.78 K, and (j) 51.55 K measured at HER and the fitting results of 
Eqs. ([2]), (EJ), and (J3J), are shown respectively. 

INS spectrum is fitted by 7 fit defined as 



rincla 



incoh / 



r\huj) = r nela (hu) + r ncoh (hu). 



(2) 



rmeia j g ^ e inelastic intensity and J mcoh is the incoherent intensity defined as follows: 
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i=l v 



(3) 



x exp 



4\n2(hu - huiif 



V7TJ-3 

r 4 
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41n2(ftu;) 5 

1 3 



(4) 



Tl+(M 2 
where q and Tj represent an intensity normalization factor and a width of the Gaussian 

or Lorentzian function, whereas hwi stands for a center of the INS peak. By fitting J fit to 

the INS spectra of \Q\ = 0.95 A -1 at various temperatures, we obtained peak positions as 

fijjji = 0.498(2) and tvuj2 = 0.607(1) meV, where hu>i and hu>2 are defined in Fig. QJb). In 

the fitting procedure, hxi\ and hu} 2 are assumed to be global parameters, and all the spectra 

at different temperatures, T = 0.71, 1.40, 2.62, 4.03, 6.05, 10.17, 30.78, and 51.55 K, are 

simultaneously fitted. The reported result [l7J expects that the excitations are observed at 

hjjj = 0.484 and 0.584 meV, which almost correspond to the present INS result. This fact 

and the Q dependences discussed later confirm that these excitations originate from the CU3 

spin cluster. 

The fitting intensity J fit for T = 0.71 K is shown by the solid lines in the Figs. [2^a) 
and[2^c). The INS spectrum is well fitted in the whole hco region. The difference between 
the experimental intensity 7 exp and the incoherent intensity / incoh is plotted in Fig. [3]^a). In 
addition, we also showed the fitting intensities 7 fit compared with the experimental results at 
T = 1.40, 2.62, 4.03, 6.05, 10.17, 30.78, and 51.55 K in Figs. E^d) -Eft). Noticeably, the INS 
spectra are also well fitted in the wide range of temperature. The background subtracted 
spectra at these temperatures are plotted in Figs. [3]^b)-l3]^h). Temperature dependences of 
the peaks become clear; the intensities of the 0.5 meV peak at 1.40 and 2.62 K are bigger 
than that at 0.71 K and then decrease when temperature increases above 2.62 K, whereas 
the intensity of the 0.6 meV peak monotonically decreases with temperature increasing. 
These temperature dependences can be explained by the Boltzmann distribution of the spin 
states as discussed later. 

The Q dependences of ftw = 0.5 meV at T = 1.80 and 10.17 K measured at HER are 
shown in Fig. HJ^a). We also depicted the Q dependences of hu = 0.58 meV at T = 1.58 
and 30.09 K in Fig. H^b). Both figures, at first glance, show only weak Q dependences nor 
temperature dependences. However, this is due to the strong incoherent contamination from 
hydrogen. To clearly see the Q dependence of the magnetic component, we subtracted the 
intensity at higher temperature from that at lower temperature. The subtracted intensities 
are shown in Figs. H^c) and@Jd). Both figures show decreasing behavior at low Q, the broad 
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FIG. 3: The experimental INS spectra measured at HER (/ cx P_/ mco11 ) anc | th e calculated intensities 
using Eq. (3) in Ref. 23 of \Q\ = 0.95 A~ x at (a) T = 0.71 K, (b) 1.40 K, (c) 2.62 K, (d) 4.03 K, 
(e) 6.05 K, (f) 10.17 K, (g) 30.78 K, and (h) 51.55 K are illustrated, respectively. 

peaks around \Q\ = 0.9 A -1 , and mostly flat behavior at high Q. They are indeed in good 
agreement with the geometry of the CU3 triangular spin cluster and the proposed interaction 
parameters between Cu 2+ ions as discussed later. 

IV. DISCUSSION 



In this section, we first determine the parameters in Eq. ([I]) from the excitations at 
hu = 0.5 and 0.6 meV corresponding to the excitations to Stotai = 3/2 from S tota i = 1/2. 
Then, we show that the low-energy excitation at huo, which is the splitting of the ground 
state quartet due to the DM interaction, is indeed observed by the high-energy-resolution 
neutron spectroscopy. Finally, we discuss the temperature dependences of the INS peaks at 
hu = 0.5 and 0.6 meV to elucidate the origin of the long life time of the spin state in the 
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FIG. 4: (a) The Q dependences of the INS peak at Huj = 0.5 meV, and T = 1.80 and 10.17 K 
measured at HER. (b) The Q dependences of the INS peak at Huj = 0.58 meV, and T = 1.58 and 
30.09 K measured at HER. (c) The temperature difference of the Q dependence at fvjj = 0.5 meV 
between 1.80 and 10.17 K and the calculated intensity, (d) The temperature difference of the Q 
dependence at fku = 0.58 meV between 1.58 and 30.09 K and the calculated intensity. 



CU3 cluster. 

The procedure to calculate the neutron scattering function was reported in Ref. 



23 



and 



23 



We also used 



we used the calculated intensity J cal (|Q|, tko) as defined in Eq. (3) in Ref. 
the magnetic form factor of Cu 2+ ions given in Ref. |24J. Since the peak widths in the fitting 
result are almost same as the resolution limited values below 10 K as shown in Fig. [7(b) 
and[Hd), and we only used the experimental data below 10 K to determine the optimum 
parameters in the spin Hamiltonian, we assume that the INS peaks have the resolution- 
limited widths in the present calculations for I cal (\Q\,hw). The temperature dependence 
will be discussed using the Boltzmann factor p^ (used in Eq. (2) in Ref. l23j) defined as 
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follows: 

ft = l^ ( 5 ) 

^ e fcB r 

j=i 
where &b is the Boltzmann constant and the labels of the energy levels are defined in 

Fig. Eft) . 

A. Parameters in the spin Hamiltonian 

To obtain the optimum parameters of the spin Hamiltonian, we performed least-squares 
fitting of J cal (|Q|, hu) to the subtracted spectra (7 cxp - J incoh ) at T = 0.71, 1.40, 2.62, 4.03, 
6.05, and 10.17 K as shown in Figs. E^a) _ EIf)- It may be noted that this time the scattering 
intensity of each peak is not an adjustable parameter, but only the overall intensity was 
optimized. The calculated scattering intensities I cal (\Q\, fko) for T = 0.71, 1.40, 2.62, 4.03, 
6.05, and 10.17 K are shown by the solid lines in Figs. EJ^a)— |3]^f); a good match to the 
observation can be readily seen. This satisfactory correspondence ensures the reliability of 
the estimated parameters. The obtained optimum parameters are as follows: 

J* 2 = J\ 2 = -4.19 ± 0.03 K, 
J* 2 = -4.67 ± 0.05 K, 

Jh = J h = J h = J !,i = - 4 - 14 ± °- 01 K > ( 6 ) 

J* 3 = J*j = -4.42 ± 0.02 K, 

D{ 2 = D z 23 = D Z Z1 = 0.66 ± 0.01 K, 

^i,2 = D i,2 = 0.55 ± 0.05 K. 

The uncertainty ranges of the obtained parameters were estimated as the standard deviation 
of the Gaussian distribution using the linear approximation. The energy levels using the 
optimum parameters are shown in Fig. QJb). The Stotai — 3/2 states are almost degenerated 
whereas the Stotai = 1/2 quartet is split into two doublets. The excitations, hoj , huii, and 
ftLU 2 ar e estimated as 0.106, 0.501, and 0.607 meV, respectively. It should be noted that the 
optimum parameters are within the 10 % difference of the reported parameters [17], which 
are Jf 2 /K = -4.49, Jf |2 /K = -4.54, Jf 3 /K = -3.91, J| 3 /K = -3.96, D z 12 /K = 0.517, and 
Df 2 /K = 0.517. In addition, the estimated excitation energies of the reported result are 
almost the same as our result. 
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FIG. 5: The experimental and calculated magnetic susceptibilities in B = 0.1 T are shown. 



The Q dependences at fko — 0.5 and 0.58 meV are also calculated using the opti- 
mum parameters; the solid lines in Figs. HJ^c) and IU(d) stand for the calculated intensities 
/ cal (|Q|, fajj). Both calculated results reproduce the temperature differences of the Q scans 
well, which confirms that the 0.5 and 0.6 meV peaks are indeed consistent with the triangular 
geometry of the spin cluster and the interactions between the Cu 2+ spins. 

To further check the reliability of the model Hamiltonian, temperature dependence of 
the magnetic susceptibility is calculated using the optimum parameters. In the calculation 
of the susceptibility, we use the reported value of the g tensor [17| in Eq. (TjQ). Figure [5] 
illustrates the comparison between observed magnetic susceptibility of the powder sample 
and the calculated susceptibility. The observed susceptibility is well reproduced by the 
calculation in the wide temperature range, again confirming the validity of the present 
parameter estimation. Therefore, the spin Hamiltonian Eq. ([T]) and its parameters Eq. (J6j) 
can explain not only the INS results but also the susceptibility measurement. 



B. Ground state splitting due to the DM interaction 

As shown in the above sections, the CU3 spin cluster most likely has energy level scheme 
summarized as Fig. QJb). However, there is a possibility that the energy levels are miss- 
assigned; i.e., the DM interaction is very small and the anisotropy of the exchange interaction 
{Jii+i 7^ Jh+i) i s l ar g e > resulting in the large splitting of the Stotai — 3/2 quartet and the 
small splitting of the Stotai = 1/2 quartet. We are now showing the definite evidence of the 
existence of the DM interaction in the model Hamiltonian by observing the excitation huo, 
or the splitting between two Stotai = 1/2 doublets. We note that the ground state splitting 
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can not be possible by the other interactions than the DM interaction. Figure [6] shows 
the INS spectra of deuterated powder CU3 [25|; S(\Q\,hui) at T = 1.5 and 30 K measured 
at DCS is integrated in the range of 0.1 < \Q\ < 1.3 A -1 . There is surely an INS peak 
at Hoo — 0.1 meV in the T = 1.5 K spectrum as expected from the HER result which is 
summarized in Fig. QJb). To4 fit the INS spectra, we use different function for I incoh (frui) in 
Eq. (|2J), instead of Eq. (jlj), as follows: 



rincoh 



(hcu) 



= c 3 - 

+ c 4 
+ c 5 



2 v / hT2 

v^r 3 
r 4 



exp 



ri + (M 2 
r 5 



41n2(W) 



ftu) 



(7) 



x 



1 — exp 



k B Tj 



The difference of the shapes of the incoherent background (Eqs. (jl]) and (J7J)) may come from 
the different resolution functions of the spectrometers. By fitting J fit of Eq. (j2j) with Eqs. (J3]) 
and (J7J) to the INS spectra, the peak position is determined as hw = 0.103(2) meV, which is 
almost exactly the same as the expected values, Jtwq = 0.102 meV of the deuterated sample 
from the DCS result and ftwo = 0.106 meV of the non-deuterated sample from the HER 
result using the excitations Tno\ and tkoz- The ground state splitting is surely confirmed by 
this high-energy-resolution experiment, and thus, we can conclude that the DM interaction 
surely exists in the Cu 3 spin cluster. The splitting of the ground state quartet is of the 
similar magnitude as that in V 3 23 ]. 



C. Spin and phonon states 



To discuss the temperature dependences of the INS peaks measured at HER in detail, the 
integrated intensities of the peaks at hv = 0.5 and 0.6 meV, C\ and Ci defined in Eq. ([3]), are 
obtained from the fitting results as shown in Fig. [2] [26]. The integrated intensities are shown 
in Figs. [7(a) andUJc), respectively. Different temperature dependences are readily seen in 
the figures. The 0.5 meV peak comes from the upper Stotai = 1/2 state to the Stotai — 3/2 
states, whereas the 0.6 meV peak comes from lower Stotai = 1/2 to Stotai = 3/2 as shown 
in Fig. [1(b). The temperature dependences of the intensities of each peaks are governed by 
the Boltzmann factor of the initial state as described in Eq. (2) in Ref. 



23l . Therefore, the 
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FIG. 6: (Color online) The INS spectra of 0.1 < \Q\ < 1.3 A" 1 at T = 1.5 and 30 K measured at 
DCS using deuterated powder CU3 are plotted. Both spectra are fitted by the same procedure in 
Fig. [2] using Eqs. ©, ©, and Q. 

peak intensity should follow the Boltzmann factor (p 3 +Pi) or (p x + ^2) as defined in Eq. (J5J) 
for the peak at ftwi = 0.5 or F1W2 = 0.6 meV as far as other perturbations are small. 

In Figs. 0a) and [7(c) , the calculated intensities, G03+P4) and (P1+P2), are compared with 
the experimental results. The observed and calculated intensities for both INS peaks are in 
agreement below 10 K, whereas they are not above 10 K. This feature is also confirmed by 
the temperature dependences of the peak widths of the INS peaks, C\ and C2 as defined in 
Eq. ([3]). Figures 0(b) and 0(d) show that the temperature dependences of the peak widths 
of the INS peaks at fko — 0.5 and 0.6 meV, respectively. Figure 0(d) shows that the width 
of the 0.6 meV peak is much broader than the resolution limit value above 10 K, suggesting 
that the relaxation time of the spin state becomes shorter above 10 K. Thus, the spin-phonon 
interaction may become relatively relevant above 10 K. It should be noted that the spin-spin 
relaxation rate I/T2 grew up under 10 K in the similar system CU3AS [121]. On the other 
hand, the calculated INS spectra for 30.78 and 51.55 K almost reproduce the experimental 
results in spite of such high temperature as shown in Figs. [3(g) and 0(h) 221] . In addition, 
the inelastic peaks were visible even at very high temperatures as 50 K. These facts imply 
that the phonons are only perturbative to the spin states even above 10 K, only giving rise 
to the broadening effect. We also note that the spin Hamiltonian Eq. (0) can reproduce the 
magnetic susceptibility up to 300 K very well as shown in Fig. Therefore, the spin state 
is only weakly influenced by the phonon (or other perturbations), and in particular, it is not 
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FIG. 7: (a) The integrated intensity of the peak at hw = 0.5 meV {c\ in Eq. ([3])) and the expected 
intensity ((^3 + p^) in Eq. ©). We use a log scale for the horizontal axis here and subsequent 
panels, (b) The width of the INS peak at tko = 0.5 meV (T\ in Eq. ([3])) and the resolution-limited 
width, (c) The integrated intensity of the 0.6 meV peak (02) and the expected intensity {p\ +P2)- 
(d) The width of the 0.6 meV peak (T2) and the resolution-limited width. 

disturbed below 10 K, at least in the neutron time scale. It is very intriguing to study the 
origin of this decoupling between the phonon and spin states, and is left for future study. 

V. CONCLUSIONS 

INS experiments have been performed on the CU3 triangular spin cluster using both 
non-deuterated and deuterated powder samples. Firstly, from the INS spectra measured at 
HER, we obtained the optimum parameters of the spin Hamiltonian listed in Eq. (|6]). These 
parameters and spin Hamiltonian can also reproduce the Q dependences of the INS peaks 
as well as the magnetic susceptibility. Secondly, we have directly observed the splitting 
of the ground state quartet due to the DM interaction at 0.1 meV; the value is almost 
exactly expected by the optimum parameters of the spin Hamiltonian determined using 
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higher energy excitations. Thirdly, the temperature dependences of the INS peaks at 0.5 and 
0.6 meV suggest that the coupling between spin and phonon in CU3 is very weak, resulting in 
the rigid spin state, or the long life time of the spin state. We, thus, conclude that they are 
the key features to explain the half-step magnetization change with the milli-second order 
hysteresis. 
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